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Abstract. Within the mode-coupling theory for the evolution of structural 
relaxation, the A4 glass-transition singularities are identified for systems of 
particles interacting with a hard-sphere repulsion complemented by different 
short-ranged potentials: Baxter's singular potential regularized by a large- 
wave-vector cutoff, a model for the Asakura-Oosawa depletion attraction, a 
triangular potential, a Yukawa attraction, and a square-well potential. The 
regular potentials yield critical packing fractions, critical Debye-Waller factors 
and critical amplitudes very close to each other. The elastic moduli and the 
particle's localization lengths for corresponding states of the Yukawa system and 
the square- well system may differ by up to 20% and 10%, respectively. 



PACS numbers: 61.20.Lc, 64.70.P, 82.70.D 

1. Introduction 

In this paper, some results shall be presented for the liquid-glass-transition diagrams 
of simple systems as obtained within the mode-coupling theory for the evolution of 
structural relaxation (MCT) [Q. This theory is based on closed non-linear equations 
of motion for the normalized autocorrelation functions for density fluctuations with 
wave- vector modulus q = \q\, 4>q(t)- Liquid states are characterized by solutions 
vanishing for large times t, 4> q (t — > 00) = 0. Glasses are characterized by solutions 
whose long-time limits are the non-vanishing Debye-Waller factors of the arrested 
amorphous structure f q = <j> q (t — > 00), < f q < 1. The f q obey the set of implicit 
equations 
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The mode-coupling functional T q reads 

E vtfAdfkU ( 1& ) 
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with coupling coefficients 

V(q, k,p) = pS q S k S p {[kc k +pc p }q/q} 2 /(2q 2 ) . (2) 

Here, S q and c q are the structure factor and the direct correlation function, 
respectively, related by the Ornstein-Zernicke equation S q = 1/ [1 — pc q ] . The structure 
factor depends on control parameters like the particle density p, the temperature T, 
and parameters specifying the dependence of the interaction potential U(r) on the 
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interparticle distance r. Let us combine the set of control parameters considered to 
a control-parameter vector V. The discussion will be restricted to parameter regions 
where S q , and hence T q , depend smoothly on V. The Debye- Waller factor f q for a 
given V is distinguished from possible other solutions of Equations for the same 
V, say f q , by the maximum property: f q > f q for all q. 

For almost all control parameters, the Debye- Waller factors f q depend smoothly 
on V. The exceptional points are referred to as glass-transition singularities V c . 
These critical points of Equations (|l]j^) are bifurcation singularities of the cuspoid 
family Ai, I — 2, 3, . . . An Ai bifurcation describes a topologically stable singularity 
that is equivalent to the bifurcation singularities of the zeros of a real polynomial 
of degree I [Q. The liquid-glass transition points are A 2 singularities located on 
smooth surfaces in parameter space. If V crosses this surface at some V c , the long- 
time limit <p q (t — > oo ) jumps from zero to the critical Debye- Waller factor f q > 0. 
There may exist surfaces of A 2 singularities within the glass. These describe iso- 
structural transitions from one amorphous state characterized by f q > to another 
state specified by a larger Debye- Waller factor f q > f q . For every A 2 singularity, a 
number A, 0.5 < A < 1, can be calculated . It is called the exponent parameter since 
it determines the various anomalous exponents entering the description of the slow 
dynamics for V near V c . The higher-order singularities, Ai, I > 3, are the endpoints 
of the transition surfaces characterized by A = 1. In the following, some properties of 
A4 singularities V* shall be discussed. Interparticle interactions shall be considered 
consisting of a hard-sphere-repulsion core of diameter d and a short-ranged attraction 
potential for r > d. The latter shall be parameterized by a typical attraction strength 
uq and a typical attraction range A. Therefore, a three-dimensional control-parameter 
space is considered: V = (ip, T,S), where tp = (ir/6)pd 3 denotes the packing fraction 
of the spheres, T = uo/(1:bT) is a dimensionless attraction strength or an inverse 
dimensionless temperature, and S = A/d is an attraction-range parameter. 

For a system of particles interacting with a hard-sphere repulsion complemented 
by a short-ranged attraction, MCT leads to two intriguing results || Q. First, the 
increase of the attraction strength T may cause a melting of the glass. This implies 
a reentry phenomenon for the transition diagram. For certain values of the packing 
fraction ip, the liquid freezes into a glass not only by cooling but also by heating. 
This prediction has been verified recently for colloidal suspensions (|, ||. In this work, 
the parameter T was varied by adding polymer to the solvent thereby increasing the 
strength of the depletion attraction. The reentry phenomenon was also established 
by molecular-dynamics-simulation studies || 0, [|, |) . Second, the existence of an A3 
singularity was predicted. The glass-glass transitions connected with this higher-order 
singularity deal with discontinuous changes of the localization mechanism from one 
caused by the repulsion-induced cage effect to one caused by the attraction-induced 
bonding of the particles. The cited work was based on Baxter's model for sticky hard 
spheres (SHS) complemented with a large cutoff wave vector q co restricting the 
wave- vector sums in the mode-coupling functional from Equation (JT^). The effective 
range of the regularized potential introduced by the cutoff can be parameterized by 
A = 7r/q co . The singular Baxter interaction leads to a large-g-tail for the direct 
correlation function c q = 0(l/q). It is this tail for 2ir/d < q < q co that causes the A3 
bifurcation. This A3 singularity disappears if g co is decreased towards 2ir/d. 

The generic scenario for the disappearance of an A3 singularity in a three 
parameter system is the existence of an A4 glass-transition singularity at some point 
V* = ((p*,T*,S*). Thus, the cited results fl], |, §, f§ lead to the following conclusion. 
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A system of particles interacting with some steep strong repulsion core complemented 
by a sufficiently strong short-ranged attraction potential exhibits an A4 singularity at 
some control-parameter point V* . The general transition diagram is organized around 
this V* . Diagrams for different models for the interparticle interaction can be mapped 
onto each other for small V — V* by a smooth invertible parameter transformation. 
In this sense, the bifurcation scenario is universal. The scenario at an A4 singularity 
has been demonstrated in some detail for the square- well system . In the following 
, these results shall be extended and compared with the ones calculated for other 
potentials that might be of interest for the description of colloidal systems. 



2. Results 



Let us specify the attraction potentials to be considered and the approximation 
theories to be used to evaluate S q . The model of SHS complemented by the mentioned 
cutoff q co shall be used with the convention T = 15/r. For the definition of the 
stickiness parameter r and the evaluation of the structure factor we follow Baxter's 
work |[o) . The hard-core Yukawa system (HCY) is used with the following convention 
for the control parameters 

U(r)/{knT) = -TeM-(r-d)/(Sd)]/(r/d), d<r, (3) 

i.e., S is chosen as the inverse of the conventional screening parameter b. The 
structure factor is evaluated analytically in the mean-spherical approximation |Q. 
Furthermore, three potentials of polynomial shape are considered: 

U(r)/(k B T) = -T[(d + Sd-r)/(dd)] {n - 1) , d<r<(l + S)d. (4) 

For n = 1, the formula describes the square- well system (SWS). The triangular- 
potential model (TRI) is obtained for n = 2. The Asakura-Oosawa system (AOS) is 
modelled by n = 3. The latter potential is obtained as small-6 limit for the depletion 
attraction acting between spheres in a dilute solvent of small polymers |T^ |. It was 
shown that the equations for the mean-spherical approximation for S q of the SWS 
can be solved analytically by an expansion in 5 |HJ . This procedure can be extended 
to treat the potentials of Equation (^) for every n. The leading order expansion 
formulae are noted in Appendix A , and they are used in the following. Let us recall, 
that for the range of length-parameters of interest, say 5 < 0.20 for the SWS, the 
next-to-leading-order result is very close to the full numerical solution for n = 1 fill . 
To determine f q for a given V, the standard iteration procedure is applied: 

/ 9 ( "V(1 - ft ] ) = ^L££ _1) ]» " = 1,2,..., /f = 1, lim„_4" ) = f q §. This is 
done after Equations (|l|,|2|) are rewritten so that the wave- vector moduli are discretized 
to M points on a grid of equal spacing h. The values h = OA/d and M — 300 have 
been used to identify the A 2 singularities on the bifurcation surface. It is cumbersome 
to identify the higher-order singularities since (I — 1) control parameters have to be 
scanned and the convergence of the iteration is slower for points at an Ai + \ singularity 
than at an Ai one. One can use the deviation of A from unity to characterize the error 
for the identification of an Ai with I > 3. In our calculations we achieved 1 — A < 10 -3 . 
At the end, the calculation at V* was repeated with h = 0.08/d and M = 1500 in 
order to check the independence of the results from the discretization approximation. 

As a representative example, figure [l] exhibits the glass-transition-singularity 
diagram of the HCY. The inner points of the lines are A2 singularities obtained as 
constant-^ cuts through the bifurcation manifold in the three-dimensional space of 
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Figure 1. Transition diagram for the hard-core- Yukawa system (HCY). The 
attraction strength T versus packing-fraction tp curves show cuts through the 
surfaces of glass-transition singularities for fixed attraction-range parameter 8 as 
indicated; cf.Equation (H). The line for 8 = 8*= 0.0173 hits the A± singularity 
(*). The line for <5 = 0.0100 exhibits a crossing point (□), and an A3 singularity 
(O)- Four further A3 singularities are marked by circles; they refer from left 
to right to 8 = 0.0167, 0.0143, 0.0125, 0.0111. The arrow points to the critical 
packing fraction ¥^35 = 0.516 for the hard-sphere system. 

control parameters V. There are three topologically different cuts. The smooth line 
shown for S = 0.0250 represents a typical cut for 6 > 5* . All points describe liquid-glass 
transitions characterized by an exponent parameter A < 1. For 5 = S* = 0.0173. . ., 
the transition line runs through the A4 singularity marked by a star. For all V ^ V*, 
the line is smooth and A < 1. For V = V*, A = 1 and the line has a continuous tangent 
but an infinite curvature. The lines for S < S* consist of two pieces as demonstrated 
for 6 = 0.0100. The part shown as a heavy line with dT c /d(p c < terminates within 
the glass in an A3 singularity marked by a circle. The second part shown as a light 
line terminates in a crossing point indicated by a square. The line of Ai singularities 
between the crossing point and the A^ singularity deals with glass-glass transitions 
and the other part with liquid-glass transitions. A2-glass-transition curves for the 
HCY for 8^5* cuts have been considered previously for the HCY f| |TJ, [l^] . 

Let us add a remark on the diagram. Obviously, for V — 0, all cuts start at the 
critical packing fraction of the hard-sphere system (HSS), (/?hss = 0.516. There is a 
characteristic range parameter Reentry so that the slope dr c /d^ c for T = is negative 
for 5 > Reentry and positive for 5 < Reentry- The latter case is exemplified in figure [j] 
by the three cuts discussed in the preceding paragraph, the former case is shown by 
the cut for S = 0.300. For 6 < <5 rcC ntry and sufficiently small T, the transition curves 
deal with melting of the glass upon increasing the attraction parameter T. Hence, for 
all liquid states with ip > ^hsS' one § e ^ s ^ ne reentry phenomenon mentioned in the 
introduction. The bonding forces create holes in the cages and this can destroy the 
particle localization for long times. In agreement with Lindemann's melting criterion, 
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the critical localization length of the HSS is about r c s = 0.0746d. If the attraction- 
potential range A is much larger than r£, i.e., if 6 exceeds a critical value <5 reentr y, the 
bonding effects cannot change the cage structure. Therefore, the reentry mechanism 
disappears for S > ^reentry- In this case, the transition curve is similar to the one 
obtained for a typical van-der-Waals system described, e.g., by a Lennard- Jones 
potential. One obtains Reentry = 0.30 and 0.117 for HCY and SWS, respectively. 
It should be noted that the transition curve shown for 5 = 0.30 terminates at the 
spinodal line defined by the divergence of S q for q = 0. Beyond this point, MCT 
equations are meaningless and calculations of density correlators cannot be performed 
there within that theory. 
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Figure 2. Attractive potentials U relative to the thermal energy k^T as function 
of the interparticle distance r relative to the hard-sphere diameter d for control 
parameters V* at the respective A4 singularity. The curves refer to the SWS ( — ), 

TRI (■ ■ •), AOS ( ), HCY ( ) and SHS ( ). The inset shows the 

cuts of the bifurcation surfaces through V* for fixed attraction-range parameter 
S = S*. 

Figure || exhibits the four regular attraction potentials defined in Equations 
for control parameters at the respective singularities. The inset shows cuts through 
the liquid-glass-transition surfaces for fixed range parameter 5* . The regular potentials 
are rather close to each other for 1.01 < r/d < 1.05. The volume of the shell 
1.00 < r/d < 1.02 is smaller than the one for the shell 1.02 < r/d < 1.04. Within 
the latter, the attraction strength decreases in the sequence SWS, TRI, AOS and 
HCY. Therefore, the critical packing fraction ip* for the onset of an iso-structural 
glass transition increases in this sequence. One finds <p* — 0.5272, 0.5307, 0.5321, 
and 0.5342, respectively. For the same reason, the maximum packing fraction of 
the liquid increases in this sequence. One gets ip m&x = 0.5293, 0.5326, 0.5340, and 
0.5367, respectively. For Baxter's model, all mode-coupling effects have been cut off 
for q > ird/8*. Therefore, the bonding effects are reduced compared to those for 
the regular potentials and the values ip* = 0.5562 and <^ max = 0.5566 for the SHS 
are larger than the corresponding values for the regular potentials. Evaluating the 
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Figure 3. Structure factor S q , critical Debye- Waller factor f%, and critical 
amplitude h q for control parameters specifying the A4 singularity V* , SWS ( — ), 

HCY ( ) and SHS ( ). The fourth curve ( ) shows the results for 

the SWS using the structure-factor up to next-to-leading order in <5 Jll| . The lines 
denoted by HSS exhibit /" and h q for the critical point of the hard-sphere system. 
The inset from the lower panel shows the variation of the exponent parameter A 
along the liquid-glass transition lines through V* from Figure ^. 

Structure factors S q at the singularities and the corresponding critical Debye- 
Waller factors /° are compared in Figure ||. Also compared are the so-called critical 
amplitudes h q . These are evaluated by a straightforward but involved procedure 
from the mode-coupling functional at the bifurcation point fil|. They quantify the 
susceptibility of the arrested structure. The increase of (f q — f q ) upon increasing 
V — V* is proportional to h q and so is the prefactor in the logarithmic decay laws 
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that are the characteristic feature of the dynamics near a higher-order singularity Jig ] . 
S q for various models mainly differ by a small shift parallel to the g-axis only. This 
shift reflects the decrease of the interparticle distance caused by the increase of <p*. 
The /° oscillate in phase with S q and the h q oscillate in opposite phase as known and 
explained for the simple HSS |Q. However, the f q are considerably larger and the h q 
are smaller than the corresponding quantities at the hard-sphere transition. This is a 
manifestation of the fact that the attractive part of the potential enforces localization, 
provided a glass state is established. Within the wave-vector region around the first 
diffraction peak, say 1 < qd < 10, the Debye- Waller factor of the SWS differs from 
the one of the HCY up to about 7% while this difference is minimal at the peak. 
The corresponding difference for the critical amplitude is about 9% which is maximal 
at the peak. The inset for the lowest panel of Figure || displays the variation of the 
exponent parameter A for the liquid-glass transitions on the cuts 6 = S*. Rescaled as 
a function of T/T*, the A cannot be distinguished on the branch T/T* < 1 dealing 
with transitions to the repulsion dominated glass. On the branch T/T* > 1, which 
deals with transitions to the attraction dominated glass, the A for the various models 
are still very close to each other. 

Figure || also exhibits results for the SWS evaluated with the structure factor in 
next-to- leading order of R eference |Tl|] in comparison with those based on the leading 



order theory explained in Appendix A. Obviously the difference is too small to be of 
interest. 

The macroscopic mechanical stiffness of liquids and glasses is quantified by the 
elastic moduli. The longitudinal modulus Ml specifies the stiffness for compressions 
and the transversal one Mr, also called shear modulus G" , the stiffness for shear 
deformations. They are defined as constants of porportionality in the linearized stress- 
strain relation. In systems with Newtonian microscopic dynamics, they determine the 
speed of longitudinal and transversal sound, respectively, via ul,t = y Ml,t / '(P m ) 
with pm denoting the mass density. For an ergodic system, the shear modulus 
vanishes, Mj, = 0. The longitudinal modulus reads M° = p{k-g,T)SQ 1 with 
So = lim 9 ^ S q . In the glass state, the moduli are larger: Ml,t = M° T + £Ml,t- 
The additional contributions are the positive long-time limits of fluctuating-force 
correlators. For the latter, MCT yields Q: 

5M L>T = p(fc B T) lim V S k S p f k f p {[kc k +pc p ]^ T } 2 (p/2q 2 ). (5) 

k+p—q 

Here are unit vectors parallel and perpendicular to q. respectively. The limit 
leads to 

roc 

« L , T = p(k B T) / dk{p[S k f k k/(27r)} 2 w h . T (k)} , (6a) 
Jo 

«*.(*) = 4 + \{kd k )c k + l -{kc' k f , (66) 
w T (k) = ^(kc' k f . (6c) 

In Figure ^ the moduli for the SWS are compared with the ones for the HCY. 
The states are chosen on the cuts 5 = S* through the respective transition surface. 
The compression modulus M° of the liquid varies smoothly throughout reflecting the 
well known increase of the compressibility with increasing attraction forces. The large 
variations of <5Ml/t reflect the strong effect of bonding potentials on restoring forces 
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Figure 4. Longitudinal elastic moduli M-^ and transversal elastic moduli Mx for 
the 8 = 5* transition line of the SWS (HCY) shown as full (dashed) lines. For 
the SWS, the structure factor from was used. In the limiting case of the HSS, 
M° = 74.1, SM T = 18.1, 5M L = 56.1. At the A 4 singularity for the SWS (HCY), 
5M T = 54.3(64.4), <5M£ = 163.4(193.6). At the maximum <5M T = 164.2(205.3), 
<5M L = 493.0(616.3). 



[Q, The same bonds are resisting shear as well as compression deformations. 

Therefore, there is no great difference in the behaviour of the two moduli. The 
contributions to SM^ due to the first two terms in Equation ( |6^) are smaller than 
the one due to the last term. Further, incidentally, these two contributions nearly 
cancel. Therefore, <5M£ differs from 3SM^. by less than 3% for T < T* and less than 
0.5% for r > r*. The universal properties of the A4 bifurcation imply the following: 
The moduli vary smoothly for all T =/= T* and at the A4 singularity there is a cubic- 
root singularity: M(T) - M(T*) oc (T - T*) 1 / 3 . The singular increase of M with T 
increasing through T* is a precursor of the discontinuous increase of M upon crossing 
the glass-glass-transition line for 5 < 5*. If T increases further, the moduli have to 
decrease since the density decreases towards zero. Therefore, the moduli exhibit a 
maximum at some value of T exceeding T* . By continuity, such maximum also occurs 
for cuts with S close but not equal to S* , as was noticed before for the shear modulus 

§0- 

The clearest quantity for demonstrating the change of the glassification 
mechanism upon increasing the attraction, is the variation of the particle's localization 
length r s as a function of T for states on the cut S = 5* (Figure [sj> . This length is 
given by the mean-squared radius of a particle's long-time probability density. The 
Fourier-transform of the latter is the Lamb-Mofibauer factor /* that is evaluated from 
an equation similar to Equation ([]]) jl| . One ends up with an equation for the inverse 
of r 2 s analogous to the equations for the moduli: 

1 f°° 

Vrl = - dk{ P S k f k f s k [c k k 2 /^ 2 }. (7) 
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Figure 5. Localization length of a tagged particle for states on the & = 8* 
transition line for the SWS (HCY) shown as full (dashed) line. At the A4 
singularity one gets r^/d = 0.0425(0.0390) as indicated by the arrows. The line 
labeled HSS marks the localization length r c Jd = 0.0746 for the HSS. 

With increasing T, r c s decreases monotonously from the value for the HSS, r c s /al = 
0.0746. The variation is smooth except for the cubic-root sigularity at T — T*. For 
strong bonding, r c s decreases to a T insensitive value of the order of the attraction 
potential range, as noted before for 5 > 5* in the HCY [Q. 

3. Discussion 

The first peculiarity of structural relaxation in systems of particles interacting with 
hard-sphere-like repulsions complemented by short-ranged attractions predicted by 
mode-coupling theory (MCT) ^ H]j namley the reentry phenomenon for the liquid- 
glass-transition diagram, has been established recently (§[ H 0, |[ This justifies to 
focus now on the more difficult task of testing the second prediction concerning the 
existence of A3 and A4 singularities. The signature of such higher-order singularities 
is the extreme stretching of relaxation curves as has been explained on the basis of 
asymptotic expansions in terms of powers of logarithms in time Jl6| . The recent 
simulation results ]?], |i| provide hints that there are such logarithmic decay laws. 
However, a more detailed analysis of the data would be required if one intends to 
arrive at compelling conclusions. The transition diagram exhibiting A3 singularities is 
organized around an A4 singularity. So far, such singularity was identified within 
a microscopic theory only for the square- well system Oj. For this system, all 
parameters and amplitudes necessary for a discussion of the logarithmic decay laws 
have been evaluated |jl8| . In the present paper, A4 singularities have been identified 
and characterized for a series of other attraction potentials in order to provide some 
information on the sensitivity of the MCT results on the microscopic details assumed 
for the interaction. 
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The critical packing fractions ip* found for the A4 singularities of the four regular 
potentials considered (Figure ||) or the related largest values for the liquid packing 
ip max exceed the critical packing fraction ip^ ss of the hard-sphere system by about 
3%. But the difference between these values for different models is about 1.4% only. 
It is not clear, whether or not experiments can discriminate between the ip* or y> max 
predicted for, say, SWS and HCY. But simulation studies for the diffusivity of the 
kind published recently [|[ ^) should be able to test our prediction ^hcy > ^sws- 

The singular potential for sticky hard spheres was introduced by Baxter jl(| 
as a certain attraction-range-to-zero limit in order to simplify equations from a 
mathematical point of view. But, within MCT, these simplifications cannot be made 
use of since the limit to zero for the range parameter 5 must not be permuted with 
the limit to infinity for the wave vectors in the mode-coupling integrals. The problem 
can be handled by introducing some wave- vector cutoff Q. This is equivalent to 
introducing a regular potential leading to a considerably larger ip* — (/?hss than found 
for the other potentials examined above. 

It is well known that the structure factor theories used in this paper lead to 
unsatisfactory results for the thermodynamic functions. Great progress has been made 
to improve on this point as is discussed, e.g., in Reference pA \. The thermodynamic 
functions are derived from the zero-wave-vector limits of correlators. However, 
different from the limit of large wave vectors, which is important to capture the 
essential physics of the glass transition, small-g effects play a minor role for the 
transition curves. As already demonstrated for the SWS in pH , artificially setting the 
structure factor input to zero for qd < 4 does not alter the transitions qualitatively. 
For r ^ r* where large q values dominate, transition lines are in accordance even 
quantitatively. 

The differences among the structure factors of the various models for wave vectors 
accessible in scattering experiments for colloids, say 2 < qd < 15, are very small if the 
states at the A± singularity are compared (Figure ||). However, the differences between 
the Debye- Waller factors and the critical amplitudes are of the order of 5%. The 
and h q determine the amplitudes of the leading-order formulae for the intermediate 
scattering functions for states near an A^ singularity in a similar manner as they 
determine the corresponding results near an A2 liquid-glass transition pq |. If they 
could be deduced from a data analysis with similar accuracy as one can do from data 
near an A2 singularity, one could discriminate, e.g., between the square- well model 
and the hard-core- Yukawa system. 

The tagged particle localization is described by the Lamb-Mofibauer factor, 
that can be measured by incoherent scattering. It is a bell-shaped function of q and 
its width can be quantified by the localization length r s : /| = 1 — (qr s ) 2 + 0(q 4 ). As 
expected from the shape of the attractive potentials at the A4 singularity (Figure |J), 
r* is smaller for the HCY than for the SWS by 10%. The MCT brings out that 
r c s becomes insensitive to T for strong attraction. For the SWS (HCY) one gets 
r c Jd = 0.022 (0.020) for T/T* > 1.8 within 1% accuracy. 

The largest difference between SWS and HCY is found at the maximum of the 
elastic moduli where the HCY is about 20% stiffer than the SWS. Taking the critical 
moduli of the HSS as reference, the maxima of SWS and HCY should be at 9 or 11 
times that value, respectively, as shown in Figure ||. 

The line of glass-glass transitions occurring for the cuts through the bifurcation 
surface for 5 < S* has some analogue in the line of iso-structural phase transitions 
from one face-centered cubic crystal to another one with a different lattice constant 
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identified for systems similar to the ones discussed above 

to what is shown in Figure [l], also the endpoint of the phase-transition line for the 
HCY exhibits a critical attraction constant T c that increases with decreasing range 
parameter S. However, the value T c for phase transitions of the SWS was found 
practically independent of 5 pq ]. Contrary to this finding for the crystal-crystal 
transition, the endpoint of the glass-glass transition line of the SWS behaves quite 
similar to what is found for the HCY. 

Summarizing, one can conclude that the MCT predictions for the dynamics near 
an A4 singularity are rather robust. At the present state of discussion, one can use 
any of the so far studied models to predict data semiquantitatively. Two reservations 
are required. First, the quality of the theories applied for the evaluation of the 
structure factors S q is not known for the high-density regime of interest. The S q 
are, however, the essential input functions for all quantitative considerations. Second, 
the range of validity of the basic MCT equations of motion are not understood, not 
even qualitatively. Quantitative results of MCT have been tested with encouraging 
outcome for the hard-sphere system |2^] , binary Lennard- Jones systems |24| , and silica 
[p5| . It remains to be tested by experiment or computer simulation whether or not 
the theory can also describe the systems studied in this paper. 
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Appendix A. Structure Factors 

For the calculation of the structure factors we use a scheme developed earlier for the 
SWS pjj . The short-ranged attraction of the potential from Equation (ji|) added to 
the hard core is treated in mean-spherical approximation and the resulting equations 
are expanded in the small parameter S. The structure factor is expressed in terms of 
Baxter's factor function Q(r) |26| , 

S- 1 =Q(q)Q(q)*, (A.l) 

p CO 

Q(q) = 1 - 2np / dr exp[iqr}Q(r) . (A.2) 
Jo 

In the leading-order approximation, one gets a shifted parabola within the hard core, 
< r < 1, 

Q(r) = V-l) + 6(r-l) + -, (A.3) 
2 n 

and a polynomial within the attraction shell, 1 < r < (1 + 5): 
K /l + 5-r\ n 

Q( r ) = _ _- , l< r <l+tf. (A.4) 



Here, the hard-core diameter d is used as unit of length. For r > (1 + 8), the factor 
function is zero. The various constants abbreviate K = T5, 

l + 2ip 12 V K -3g , Gif K 



(1-^)2 (l-<p) n ' 2(1-^)2 (i-p) 
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The limit 
Reference 



can be carried out in Equations (A.3-A.5) in the sense of 



It yields the result for S„ used in 



However, this procedure 



appears unjustified. Within the mode-coupling functional, Equation (JT^), one needs 
values of S q also for large wave- vectors. But the limits q — > oo and 5 — > cannot 
be permuted. Introduction of the cutoff is one way to redefine the model so that the 
MCT remains meaningful. A more appealing procedure is to consider the 5 expansion 
from Reference JlT[ and define as Baxter limit the one where terms of order S 1 are 
neglected c ompared to those of order 6°. But this is just the approximation defined by 
Equations ( A. 3 - A. 5 ) that was used in this paper. It is in no respect more complicated 
to handle than the approximation proposed originally fnj . 
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